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Abstract. The Einstein evolution equations are studied in a gauge given 
by a combination of the constant mean curvature and spatial harmonic 
coordinate conditions. This leads to a coupled quasilinear elliptic-hyperbolic 
system of evolution equations. We prove that the Cauchy problem is locally 
[ strongly well-posed and that a continuation principle holds. 

For initial data satisfying the Einstein constraint and gauge conditions, 
the solutions to the elliptic-hyperbolic system defined by the gauge fixed 
Einstein evolution equations are shown to give vacuum spacetimes. 
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^ '. 1. Introduction 

In order to construct solutions to classical field equations with constraints, 
such as the Yang-Mills and Einstein equations, it is often necessary to rewrite 
the system, either by extracting a hyperbolic system, or by performing a gauge 
fixing. The gauge fixing may result in a hyperbolic system, as is the case for 
example using Lorentz gauge for the Yang-Mills equations, or space-time har- 
monic coordinates for the Einstein equations. For discussions of hyperbolicity 
O . and gauge choices for the Einstein equations, see [11, 6]. 

On the other hand, there are interesting gauge choices which lead to a coupled 
elliptic-hyperbolic system, such as the Coulomb gauge for Yang-Mills, which 
was used in the global existence proof of Klainerman and Machedon [10]. For 
the Einstein equations, the constant mean curvature gauge leads to an elliptic 
^ ■ equation for the Lapse function, and to an elliptic-hyperbolic system for the 

^ I second fundamental form kij, cf. [4]. 

In this paper we introduce and study a gauge condition for the Einstein 
equations, which is a combination of constant mean curvature gauge and a 
spatial harmonic coordinate condition. This leads to an elliptic-hyperbolic 
system where the hyperbolic part is a modified version of the Einstein evolution 
equations, and where the elliptic part consists of the defining equations for Lapse 
and Shift. 

1.1. The gauge fixed vacuum Einstein evolution equations. Let M be 

a compact, connected, orientable C°° manifold of dimension n > 2 and let 
M = M X M. We define t : M — > M by projection on the first component. We 
will consider Lorentz metrics g on Af so that the level sets of t, Alt = {t} x M 
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are Cauchy surfaces. When there is no room for confusion we write simply M 

instead of Mf. 

Given (M,g), let T be a time— like normal to Mt, let the Lapse function N 
and Shift vectorfield X be defined by 

dt = NT + X, 

and assume AT > so that T is future directed. 

Let {ei}f^^, be a time independent frame on M and let {e*}^^]^ be its dual 
frame. Let eo = T and let {ea}a=o = {eo, ei, . . . , e„} so that {ea}a=o is an 
adapted frame on M (which is neither global nor time-independent), with dual 
frame {e"}2=o- some cases we will use local coordinates x*, and use the 
notation di = d/dx^, for the coordinate frame and the corresponding first de- 
rivative operators. In the following we will use frame indices, unless otherwise 
specified, and let greek indices take values in 0, 1, . . . ,n while lower case latin 
indices take values in 1, . . . , n. 

The Lorentz metric g is of the form 

g = -N'^dt (g)dt + gij{e' + X'dt) {e^ + X^dt), (LI) 

where g = gije^ is the induced metric on M. The second fundamental form 
A; of M in M is given by 

% = -\^Tgij = - ^N~^{dtgij - CxQij), 
where C denotes the Lie-derivative operator. The vacuum Einstein equations 



Rap = (L2) 

can be written as a system of evolution and constraint equations for (5, k). The 
vacuum Einstein evolution equations are 

dtQij = -2Nkij + Cxgij, (L3a) 

dthj = -ViVjN + N{Rij + tikkij - 2kimk"]) + Cxhj. (L3b) 

and the vacuum constraint equations are 

R-\k\^ + {txkf = (L4a) 

VjtrA; - = 0. (L4b) 



A solution to the Einstein evolution and constraint equations is a curve t 1— > 
{g,k, N, X) which satisfies (1.3,1.4). Assuming sufficient regularity, the space- 
time metric g given in terms of (g,N,X) by (1.1) solves the vacuum Ein- 
stein equations (1.2) if and only if the corresponding curve {g, k, N, X) solves 
(1.3,1.4). The system (1.3,1.4) is not hyperbolic, and to get a well-posed evolu- 
tion problem we must modify the system. We will do this by fixing the gauge. 

Let 5 be a fixed C°° Riemann metric on M with Levi-Civita covariant de- 
rivative V and Christoffel symbol Ff^-. Define the vector field V'' by 

V'' = g'^e''{Viej-Viej), (1.5a) 

or in local coordinates, 

y^ = 5^^(r^.-ff..). (1.5b) 
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Then —V'' is the tension field of tlic identity map Id : {M,g) {M,g), so that 
Id is harmonic exactly when V'^ = 0, see [5] for background on harmonic maps. 

The constant mean curvature and spatial harmonic coordinates (CMCSH) 
gauge condition is given by the equations 

tr^A; = t (constant mean curvature) (1.6a) 

= (spatial harmonic coordinates) (1.6b) 

Remark 1.1. In this paper we will restrict our attention to the homogeneous 
gauge conditions given above. However, it seems likely that the ideas presented 
here can be generalized to include gauge source functions, with the gauge con- 
ditions (1-6) replaced by for example 

tTgk = t + f°, (1.7a) 
V'' = /^ (1.7b) 

where f^jf'^ are a function and a space-like vector field on M, independent of 
the data. 

Let the second order operator Ag be defined on symmetric 2-tensors by 

^ghij = —Vm{g"'''figVnhij), (1.8) 
Ms 

where fig = \/det g is the volume element on {M,g). Using the identity 
^ m{g"^^ IJ-g) = ~^^IJ-gi ^g^ij may be written in the form 

A. u.. — V V h — V"^V h ■ ■ 

In particular, if the gauge condition F = is satisfied, Aghij = g'^^Vm^nhij- 
A computation, cf. section 3, shows that 

Rij = - ^^g9ij + Sij[g, dg] + 5ij. 
where the symmetric tensor Sij is defined by 

Sij = l{^^Vj+VjVi), (1.9) 

and Sij[g,dg] is at most of quadratic order in the first derivatives of gij. Thus 
the system gij i-^ Rij — 6ij is quasilinear elliptic. 

In order to construct solutions to the Cauchy problem for the system consist- 
ing of the Einstein evolution and constraint equations (1.3,1.4) together with 
the gauge conditions (1.6), we will consider the following modified form of the 
Einstein evolution equations, 

dtgij = -2Nkij + JCxgij, (1.10a) 
dtkij = -ViVjN + N{Rij + tvkkij - 2kimk'] - Sij) + Mij, (1.10b) 



4 



L. ANDERSSON AND V. MONCRIEF 



coupled to the elliptic defining equations for N,X, needed to preserve the im- 
posed gauge conditions, 

-AN + \k\'^N = 1, (1.11a) 

AX' + R'fXf - CxV' = + 2V"^X")e^(V,ne„ - V^en) 

+ 2V'^Nki^ - V'Nk^. (1.11b) 

If Sij = 0, in particular if F'^ = 0, then (1.10) coincides with the Einstein 
vacuum evolution equations (1.3). The vacuum Einstein evolution equations 
in CMCSH gauge is the coupled system (1.10-1.11). In view of the fact that 
gij I— > Rij — 5ij is elliptic, the system (1.10) is hyperbolic, and the coupled 
system (1.10-1.11) is elliptic-hyperbolic. 

2. The Cauchy problem for quasi-linear hyperbolic systems 

In this section we prove that the Cauchy problem for a class of quasi-linear 
hyperbolic evolution equations, which includes coupled elliptic-hyperbolic sys- 
tems of the form (1.10-1.11), is strongly well-posed. The techniques used are 
not new, cf. [12, 17, 3, 2] for treatments of various aspects of the problem for 
classical quasi-linear hyperbolic systems. The methods of Kato [8] for general 
quasi-linear evolution equations can presumably be used to prove the results 
stated here. However, in view of the abstract nature of the techniques involved 
in the approach of Kato, we have decided to give a reasonably complete treate- 
ment of the Cauchy problem for the class of equations that is of interest from 
the point of view of applications in this paper and its sequel [1]. 

Let A[g] be the ellipticity constant of g, defined as the least A > 1 so that 

A-^g{Y,Y)<g{Y,Y)<Ag{Y,Y), ^YeTM. (2.1) 

In the following, all norms and function spaces will be defined with respect to 
g. Let g be defined in terms of g, N, X by (1.1). Let 

A[g] = A[g] + \\N\\loo + \\N-^\\l^ + \\X\\loo. (2.2) 

For a curve of metrics t g{t), it is convenient to define 

A(T) = sup A[g{t)]- (2.3) 
te[o,T] 

We write D for first order spatial derivatives. The action of D on tensors 
is defined using the covariant derivative V w.r.t. g. Let Dg be the first order 
spatial derivatives of g. Then 

\Dg\ + \DN\ + \DX\<C{A[g])\Dg\, 

\Dg\<C{A[g])i\Dg\ + \DN\ + \DX\), 

where | • | denotes the pointwise norm. 

We recall some definitions and facts from analysis which will be needed in 
the proof of local existence. Let = g^^ViVj be the Laplace operator defined 
with respect to the background metric g, acting on functions or tensors, and 
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let {D) = (1 - Ag)i/2. Let W^'P denote the Sobolev spaces and let = W^'^. 
Then W'P = {DyLP, for s € M, 1 < p < oo, with norm 

IIuIIh/^.p = ||(D)*u||lp. 

In case s is a non-ncgativc integer, W'P, 1 < p < oo is the closure of C°°{M) 
w.r.t. the equivalent norm Y2\k\<s II^'^^IIlp- We will without further notice use 
the same notation for spaces of tensor fields as for spaces of functions on M. 

For / C M an interval, we use the notation F{I; W^'P) for the space of 
curves of class F with values in W^'P. Spaces which will be used are F = 
C, C°'^, C-', L°°, L^, W^'^, where C°'^ denotes the space of continuous functions 
with one (time) derivative in L°°. 

We use the notation OV^ for pseudo-differential operators with symbol in 
the Hormander class Sf ^, see Taylor [17] for details. For P G OP^, s eR, 

\\Pu\\wr,p < C\\u\\wr+s,p, for r e M. (2.4) 

In particular, {Dy G OP^ and D G OP^. 

The following basic inequalities will be used. Assume 1 < p < oo. 

1. Product estimate I (Kato and Ponce [9, Lemma X4], [17, (3.1.59)]). If 
s > 0, W^'P n is an algebra, and the inequality 

||nf||H/s,p < C(| |u| I Ivi/'s.p + ||ti||vK«'p||'f ||l°°) (2-5) 

holds. In particular if s > n/p, then < C||M||vF«.j'||i'||iy''.p- 

2. Product estimate II (special case of [13, Theorem 9.5 3], see also [17, §3.5]). 
Assume tj > 0, i = 1,2, some ti > 0. Then for s < min{ti,t2,ti+t2 — n/p), 
(where the inequality must be strict if sonic ti = n/p), 

\\uv\\wp < C\\u\\^rt^,p\\v\\^rt2:p. (2.6) 

3. Composition estimate ([17, §3.1], see also [15, Theorem 1]) Let 1 < s < fi, 
and suppose F G C^(M) with F(0) = 0, where denotes the Holder 
space. Then for u G D 

||F(u)||vi/s>p < C||n||vi/s>p(l + lliill^^^). 

4. Commutator estimate I ([17, Prop. 3. 6. A]) Assume P G OP'^, s > 0, 
a > 0, then 

II [P, n]t;||vi/<^,p < C(||L'u||i;,^||w||^^s-i+ff,p + ||n||vi/s+<T,p||t'||L°°)- (2.7) 

5. Commutator estimate II ([17, (3.6.2)]) Assume P G OP^. Then 

||[P, u]t;||LP < C||DM||L<x>||f IIlp- 

Restricting to Sobolev spaces of integer order, the above inequalities can be 
proved using the classical methods of calculus. 

We next introduce the class of nonlinear evolution equations which will be 
considered. As our application is to the Einstein evolution equations, we will 
consider symmetric 2-tensors on M as the unknowns, but it should be stressed, 
that the proof generalizes essentially without change to sections of general vector 
bundles over M. 
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We will think of the symmetric 2— tensors Uij , Vij as sections of the vector 
bundle Q of symmetric 2-tensors over {M,g) with fiber inner product {u,v) 
given by 

{u, v) = UijVkig''^ g^^ 

and the corresponding norm \u\ defined by |u| = (u, -u)^/^. The natural fiber 
inner product on derivatives is 

{Vu,Vv)g = {VmU,VnV)g'^'', 

with corresponding norm iV^lg. The covariant derivative is metric, 

Y{u, v) = {Vyu, v) + {u, Vy -y) , 
where (Vyu)ij = Y"^VmUij, and the rough Laplacian on Q, defined by 

cf. (1.8), is self-adjoint with respect to the natural inner product on Q, 

{VU,Vv)gfJ,g = - {AgU,v)lIg. 

'M Jm 

The curvature on Q can be computed in terms of the Riemann tensor R'jf^i of 

Let the operator L with coefficients given by {g, N, X) be defined by 

L[g, N, X] = f,^*^ - - ^ V (2.8) 

We will use the notation U = {u,v), W = x H'-\ W^'°° = x L~. 

Write L[U] for L given by (2.8), with {g,N,X) = {g,N,X)[U], and consider 
Cauchy problems of the form 

L[U]U = r[U], U\^^^ = ZY°. (2.9) 

Define the space C^{W), 1 < < [sj, to be 

4(H^) = no<,<fc-iC^aO,T];H^-^y 

Definition 2.1. A number T > is called a time of existence in for 
the Cauchy problem (2.9) if there is a unique solution U G C([0,T];7^*) n 
C^([0,r];H^-^) to (2.9). The maximal time of existence in W for 2.9 is 

T+ = sup{r :T is a time of existence for (2.9)}. 

The Cauchy problem (2.9) is called strongly locally well posed inC^iTi^) if 

the solution map —>-U is continuous as a map 

for a time of existence T = T{U^) > 0, which depends continuously onU^ G Ti^ . 
The continuity ofU^^U is called Cauchy stability. 

The following definition states the regularity properties of L,J^ which will 
imply that the Cauchy problem (2.9) is strongly locally well-posed. 
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Definition 2.2. Let V C be connected and open. The system L\l/{]U = ^[W] 
is called quasi— linear hyperbolic in V, if the maps 

U ^ h[U] = {gp\,N[UlXp\), V ^ H' (2.10a) 

U ^ J^[U], V^H' (2.10b) 

are defined and continuous, g = g\U] defined in terms of {g, N, X)p{] satisfies 
A{g) < oo for U G V, and there is a continuous function Cl = Cl{U^) such 
that for each the following holds. 

1. By^^(U^) C V, where By^^{U^) is the ball in W of radius 1/Cl, cen- 
tered at . 

2. The maps (2.10) are Lipschitz on By^^lU^) with Lipschitz constant Cl 
w.r.t. Tf , 1 <r < s, explicitly 

\\h[Ui]-h[U2]\\Hr <CL\\Ui-U2\\-Hr 

Wi]-rpl2]\\nr < CL\\Ui-U2\\nr 

for allUiM2 G B{^^^{W). 

3. The maps h, T have Frechet derivatives Dh, DJ^ satisfying 

\\Dh[U]U'\\Hs-i < Cl\\U'\\^s-i 
\\DJ^[U]U'\\ns-i < CL\\U'\\ns-i 

for allUe Bli^^{U^). 

4. Let m be an integer, 1 < m < [sj — 1. For all integers j, 1 < j < m, 
{Qi=i, > 1; Ei=i^i = ihe Frechet derivatives D^h[U], D^J='[U] of 
order j are Lipschitz functions from By ^^(U^) to the spaces of multilinear 
maps 

j 

i=l 
j 

i=l 

respectively. We call m the order of regularity of (2.9). 

Remark 2.1. The order of regularity determines the regularity of the solution 
w.r.t. time. For strong local well-posedness, it is sufficient to have order of 
regularity m = 1. 

The following is the main result of this section. 

Theorem 2.3. Let s > n/2 + 1. Assume that (2.9) is quasi-linear hyperbolic, 
regular of order 1 < m < [sJ — 1, m V C H^. Then the following holds. 

1. (Strong local well-posedness) The Cauchy problem (2.9) is strongly locally 
well-posed in C^~^^{7{^) with time of existence which can be chosen de- 
pending only on M, A[g[U%\\Dg[U%\L-- , \\9p(%\H-, and the 
constant Cl = Cl{U^) in Definition 2.2. 
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2. (Continuation) Let be a maximal time of existence for (2.9). Then 
either the solution leaves V at or = oo or 

limsupmax(A[5], ||i:)^||Loo, ||5t5||Loo,CL) = oo. (2.11) 
t/'T+ 

The rest of this section is devoted to the proof of Theorem 2.3. In subsec- 
tion 2.1, the basic energy estimate, Lemma 2.7 is proved. Local existence and 
uniqueness is proved in subsection 2.2, Cauchy stabiUty is proved in subsection 
2.3 and the continuation principle, point 2 is proved in subsection 2.4. 

2.1. Energy estimates. In this subsection, fix T > 0, s > 1, let M = [0, T] x 
M, and assume that g G L°°([0, T]; n C°'i([0,T];L°°) and N,N-\X e 

L°°([0, T]; Unless otherwise stated, all constants in this subsection de- 

pend only on T and A(r). For the applications in this paper, there is no loss 
of generality in assuming, in this subsection, that all fields are C°° on M. We 
will consider the linear system 

L[g,N,X]U = J^, Ul^^=U\ (2.12) 

where L is given by (2.8) and 



V 



Given g,X, let p be defined by 

P = -\{dtg-Cxg), (2.13) 

so that N~^p is the second fundamental form of M in {M,g). Then under the 

present assumptions, p G L°°{M). 
Define the energy £ = £{t,U) by 

^ JMt 

Lemma 2.4. Assume that U G L°°([0, T]; p C°'^([0, T]; 7^^) is a solution 
to (2.12). Then 

\dt£\ < C{£^/^mn^ + (1 + \\p\\l^)£), (2.14) 

where C = C{A\g\). 

Proof. A computation shows 



dt£ = [ ({u, Nv + F„) + {Vmu, VnFu)g'^' 

JMt ^ 

+ {VmU,nrnU)X^g"'^ + {v,F^)') Pg 



Mt 



1 



+ I [ {ViU,Vju)p'^ - -{\u\' + \Vu\'g + \v\')tTp ) pg, 



where the curvature term TZmU is given by 

T^rn'U'ij — +-R irn^lj ~^ ^ jrn^il' 

An application of the Schwartz inequality gives the result. □ 
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Using (2.14) to estimate dtE^I"^ and integrating the resulting inequality gives 

_^V2(0)| < C (\\nL^ilO,T];n^)+l\l + \\pmL-)£'/\t)dt^ • 

(2.15) 

An application of the Gronwall inequality gives 

£'/\T) < (7e^/o^llpWlk-'ii(^V2(o) + ||J^|Ui([o,t];W1))- (2-16) 
Higher order regularity is proved by estimating Es{T) = Es{T;U) defined by 
E,{T) = \\U\\L^^lo,T];n^), s>l. (2.17) 

We have 

Es{T;U) < CEi{T; {Dy'^U), s > 1. (2.18) 
The inequality (2.16) gives 

Lemma 2.5. AssumeU G L~([0, T]; 7^2)nC°'i([0, T]; H^) solves (2.12). Then 
with p given by (2.13), 

Ei{T-U) < Ce^^oMmL^dt^^^^Q.^^ + ||^||l1([o,t];H1))- (2-19) 

□ 

In order to derive higher order energy estimates, we will apply (2.19) to 
{DY~^U using the identity 

L{DY-^u = [L, {DY-^]u + {Dy-^j^. 

In order to do this we must estimate the commutator [L, {DY~^Y^. This is 
done in the following Lemma. 

Lemma 2.6. Let s > I. Assume {g,N,X) e n W^'"^ , U e W n W^'°° . 
Then 

\\[L, {Dr-']U\\n^ < C{A[g]){M\w^>^\mw + ||^||//H|W|l>Vi.-)- (2-20) 
Proof. By construction, [dt, (D)] = 0, and hence 

In the following we will write D for a first order operator with smooth coeffi- 
cients, such as given by for example V G OV^, so that Vx = ED with B of 
order zero. Recall that [D, {D)^^^] G CP*"^. In case s is an integer, then V*^^ 
may be used instead of {D)^~^ in this proof. We need to estimate the following 
quantities: 



{Dr-\N]v\\Hr, \\[{Dr-\Vx]u\ 



\\[{Dr-\NAg]u\\L2, \\[{Dr-\Vx]v\\L2. 

We will treat each term separately. For the first term, the commutator estimate 
gives 

\\[{Dy-\N]v\\Hi < C\\DN\\loo\\v\\hs-i + ||Ar||^.||^;||ioo. 
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The identity 

[BD, {Dy-'^]u = B[D, {Dy-^]u + [B, {Dy-^]Du, (2.21) 
gives using the product and commutator estimates, 

\\[{Dr-\Vx]u\\Hi < C{\\X\\w,,^\\u\\hs + \\X\\hs\\Du\\l^). 
This takes care of the second term. Next, the identity 

[{Dy-\BD]v = {Dy-^B, D]v + [{Dy-\D]Bv + [D{Dy-\B]v 

gives 

\\[{Dy-\VxML^ < C7(||X||vKi,oc|H|i/.-i + \\X\\hs\\v\\loo), 

which takes care of the fourth term. Finally, for the third term, we may in view 
of (1.8) write the second order operator P = NAg in the form 

Pu = AD'^u + BDu. 

With D = V, we have 

AD'^U = Ng'^'^VrrNnU, 

BDu = -NV'Viu = -Ng^-{r'^^ - t'^JViu, 

with AeH^n W^'°°, B G H^-^ n L°°. 

The term [{D)'^~^ , BD]u is can be estimated in by expanding the com- 
mutator and using the product estimates to get 

\\[{Dy-\BD]u\\L^ < C{\\B\\l'^\\u\\h^ + \\B\\hs-i\\Du\\l'^). 

The identity 

[AD'^, V = A[D\ {Dy~^]u + [A, {DY-^D]Du + {{DY~^ D A){Du) , 

(2.22) 

together with the commutator and product estimates gives 

\\[AD\ {DY-M\l^ < C{\\A\\wi,oo\\u\\hs + \\A\\hA\Du\\loo). 
We now have 

||[P,(L>)^-iHU2 <C((||A||h.i,oo + ||B|Uoc)|H|^^. 

+{\\A\\hs + \\B\\hs-i)\\Du\\l^.) 

This gives 

||[(D)^-i,iVA,]u||i. <C7(A[5])((|MlH^i.oo + ||7V||H/i,-)||n||j^. 

+{\\g\\H^ + \\N\\hs)\\Du\\l^). 

Collecting the above and using A[^], ||^||, H-D^H for the terms involving norms 
of g, N, X, gives the result. □ 

Using the Gronwall inequality. Lemmas 2.4, 2.5 and 2.6, gives the following 
higher order energy estimate for the linear system (2.12). 
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Lemma 2.7. Assume g £ L°°([0, T]; i?''nTy^'°°), dtg G L°^i[0,T]; H''-'^nL°°), 
r > 2. LetU = {u,v) he a solution to (2.12) satisfying u G L°°([0, T]; iJ'' n 
W^'°°)nC^'^{[0,T];H''-^), V G L'^{[0,T]-H'-^ nL'=^)nC^'^{[0,T];H'^-^), r > 
2. Then, for 1 < s < r, there is a constant C = C{T, A{T)) so that 

+ f \\g\\Hs\Mw^,o.dt). (2.23a) 
Jo 

In particular, if s > n/2 + 1, 

Es{T-U) < Ce^Io^M\H^-^+MMd\EMl^) + \\HLH[o,T];n^})- (2.23b) 

□ 

To prove (2.23b), note that if s > n/2 + 1, the quantities ||p||/^oo, ||D^||ioo, 
and \ \U\\ yy l,oo , cir6 dominated by II^H^s, and Es{T;U). 

2.2. Local existence and uniqueness. This subsection is devoted to the 
proof of local existence and uniqueness for (2.9). The proof is an iteration 
argument, following [12]. 

Note that by Definition 2.2, Cl = Cl{U^) is continuous in U'^. Further, 
A[5[W]], ||Li5[W]||i,oo can be estimated in terms oi A[g[U%, {{DglU^lL-o and Cl- 
In the rest of this subsection, unless otherwise stated, constants will depend 
only on M, A[g[U%\\Dg[U°]\\Loo , \\U°\\w, MU^WIh^, and the constant Cl = 

CLm. 

Let R = 1/Cl- Let {Z^^}5^=o C C°° fl B^^^(U'^) be a sequence of approxima- 
tions of W° given by smoothing. We will construct a sequence of approximate 
solutions 

for some T^, > 0, to be chosen, with initial data for Um given by U^. 

Set Uo{t) =Uq,J^o = 0, and let Lq = L[gQ, A^O; ^o] be the time-frozen version 
of L defined by setting {goit), No{t), Xo{t)) = {g[U^],N[U^],X[U^]). 

For m > 0, let 

{gm,Nm,Xm) = {g[Um], N[Um], X[Um]), 

and for m > define Um+i to be the solution of the linear Cauchy problem 

LnMm+l = Fm, Um+l\^^Q = Z^m+1- (2.24) 

The existence of solutions for the Cauchy problem for linear hyperbolic problems 
with smooth coefficients is standard, the proof given in [16, Theorem 3.3] is 
easily adapted to the present situation. Suppose {l^m'}m'=o ^ sequence of 
solutions to (2.24) with Um' taking values in B^^{U^). By point 3 of Definition 
2.2, we get 

\\dtg[Umm\Hs-^ < CL||5tW„,||^.-i. (2.25) 
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Let pm = — \{dtgm — Cxm9m)- From (2.24) it follows, using the energy estimate, 
that llSfZ/Zmllw-i is bounded by a constant depending on Cl and hence in view 
of (2.25) we get a bound on Hpmlli?"-! + Hymll//"- Given this estimate, it follows 
from Lemma 2.7, that there is a constant Cr < oo, so that as long as {Um'}m=o 
takes values in Bj^{U'^), the energy estimate 

Es{T;Um+l) < CR{Es{0;Um+l) + \\^m\\L^{[0,T];W)) (2.26) 

holds for Um+i, for T < 1. The restriction T < 1 is made so that Cr does not 
depend on T. We can choose the sequence {W^j^^Q such, that 

W° G B%/^{U°) for m > (2.27) 
CR\\Um{0,-)-Um'{0,-)\\ns < R/4 for m,m' > 0. (2.28) 

We will prove, for a < 1 to be chosen, convergence for this sequence in 
L°^{[0,T4;H')nC^^\[0,n];n'-^) to a limit 

by the following steps: 

1. {Um}^^o C i°°([0,r,];S^(WO)), by Lemma 2.8 

2. {Z^m}m=o is a Cauchy sequence in L~([0, T*]; Ti^) n C^'\[0,n];n^), by 
Lemma 2.9, with limit U satisfying 

by Lemma 2.10. See (2.30) below for the definition of C^,. 

3. W is a solution to (2.9). 

4. U e C{[0,T4;H') nC\[0,T4;H'-^), Lemma 2.11. 

Lemma 2.8 (Boundedness in high norm). There is a time T* > such that 
{W^}CL~([0,T,];5^(W0)). 

Proof. We must prove that there is a > so that if ^4^ G L°°{[0, T*]; S^(Wo)), 
then Um+i G L~([0, T*]; 5^(Z^°)). In order to do this we consider 

Lm{^m+l — Uq) = — LrrMo 

It follows from (2.26) that 

\\i^m+l -^o\\l'^{[0,T];W) < Cr (||W^+i -Uo\\w + 

\\^m\\L^{[0,T];n^) + \\LrrMo\\L^{[0,T];H')) 

and we see using (2.28) and the fact that ||/||Li([o,t]) ^ ^ll/llL°°([o,t])) that there 
is a < 1 so that if G L°°{[0,T^]; B%(U^)), then 

\\^m+l - Wo||L°o([0,r*];H'') < (2.29) 

By construction Uo{t) = Uq and hence it follows from (2.27) that Um+i G 
L°°([0, T*]; B^(W°)). This completes the proof of Lemma 2.8. □ 

Lemma 2.9 (Convergence in low norm). There is a time T* > 0, so that {Um} 
is Cauchy in L°°([0,r*];7^i) n C°'i([0,r*];H°). 
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Remark 2.2. In order to handle the cases with s < 3, we show in Lemma 2.9 
that {Um} is Cauchy in L'^{[0,T^];Ti.^)r\C^'^{[0,T^];Ti.^), in particular we have 
control of dtVm only in . In case s > 3, this can he avoided and'H},Tf' can 
in the rest of this section he replaced hy Ti?,Ti}. In particular if s >3, dfVm is 
Cauchy in L°°{[0,T^];L^). 

Proof. Let be as in Lemma 2.8. Let T <T^. We compute 
and hence 

ll^m+l - Z^m'+l||L°°([0,T];Wi) < Cr (||Z^m+l -^m'+lWn^ 

+ W^m - •^m'||Ll{[0,T];Wl) 
+ ||(-^^m - -^m')^m'+l||Ll([0,T];Hl)) 

Using the Lipschitz property of the map U {g, N, X, J-) we see that by 
possibly decreasing T^, we get for all m,m' > mo, 

\\l^m+l -^m'+l||L°°([0,r.];Hi) < Cr\\U^_^_i - U^'+i)\\n^ 

+ l^^l^ra — Wm'||L°°([0,T.];Hi) 

As is Cauchy in H} , we can by thinning out the sequence \U.^ get 

m 

Let = Cr||^^+i -U^W-ni, m > and am = \\Km - ^m-i||L°°([o,r.];Wi), 
m > 1. Then, 

Om+l < Pm + ^Clm, m>l. 

Solving the difference equation 
gives 

oo oo 

^ = 2 ^ /3„ + 2ai 

m=l m=l 

and which shows {Um} is Cauchy in L^{[0,T^];7i^). 

In the equation of motion, which determines dfUm, there will occur terms 
of the form g^-'didjUm and X^djVm- In order to show that dtUm is Cauchy in 
L°°{[0,T:t:];H^), we need to show that these terms are Cauchy in II~^. To see 
this, recall that multiplication is continuous x ^ for s > n/2+ 1, cf. 
product estimate II. In view of the fact that II~^ is the dual to H^, this implies 
that multiplication defines a continuous map x — > for s > n/2 + 1. 
It is now clear from the mapping properties of {g, N, X, J^) and the equation of 
motion, that {Um} is Cauchy in L~([0, T^];H^) n C^'^QO, X,]; □ 
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The dual space to -ff* is Let {(p,u)-s,s denote the duahty pairing of if ^ 

with Then \{(j),u)^sj < MIh-AMh- ■ Define the space Cu,{[0,T]; H^) of 
weakly continuous functions on [0,T] with values in H^, i.e. 

CU[0,T];H') = {u:{cl),u)^s,s(^C{[0,T]), for ah G i/-*} (2.30) 

Lemma 2.10 (Weak convergence). Let {umj^^i C C{[0,T]; H""), s > be a 
bounded sequence and assume {um} is Cauchy m L'^{[0,T];H' ) for some s , 
< s' < s. Then there is a u e Cyj{[0,T]; H^) D L'^i[0,T];H^) so that for all 
4> G , {(f), Um)-s,s {(f>j '")-s,s uniformly in L°°[0, T] as m oo. 

Proof. Let 4> G be arbitrary and fix e > 0. Let C be a constant so that 
||'Um.||_f/s < C for all m. Let u be the limit of {um} in L°°{[0,T]; ). By the 
uniform bound on ||'Um||L°°([o,r];ii''')) we find u G L°°([0, T]; i?*). 

Recall H~^' is dense in H~^, hence there is a, (f)' E H~^' so that ||i?!> — 
4>'\\h-' < e/3C. Since tim u m , wc may choose m large enough so 
that ||(/>'||^_s/||uot — ■u||^s/ < e/3. Then at t G [0,?"], 

\{<l),Um)-s,s- {(l>,u)-s,s\ < \{4>- (f>',Um)-s,s\ 

+ \{(l>',Um - u)-s',s'\ + \{4>' - (t>,u)-s,s\ < e 

As e was arbitrary, we find that the sequence {(p,u„i)-s,s G C{[0,T]) converges 
uniformly to {(p,u)-s,s e C{[0,T]). Hence u G Ci„([0, T]; i?*) as ^ G H'' was 
arbitrary. □ 

By Lemma 2.10, Um-^U e C^([0, T^\,W). Next we prove that U solves (2.9). 
It is clear from the construction that U{Q) = U^, so we need to show that 
L[U]U — J^[U] = 0. Consider the sequence {Um} defined by (2.24). We compute 

L[U]U - T[U] = {L[U] - L[Um])U + L[Um]{U - Um+i) + {:F[Um] - J'iU]). 

By Lemma 2.9, {Um} converges to U in L°°([0, T*]; 7Y^) n C°'i([0, T*]; W°). By 
the Lipschitz property of U ^ {g,N,X,J^), the first and third terms in the 
right hand side tends to zero in L°°{[0,T:^];H^) as m — > oo. It follows from the 
definition of L and the discussion in the proof of Lemma 2.9 that the second 
term tends to zero in L°°([0, T*]; TiP). However, the left hand side is independent 
of m and hence equals zero. This proves that ^ is a solution to (2.9). 

At this stage we know that t (U,dtU) is weakly continuous and that U is 
a solution to (2.9). In order to prove U G C{[0,T^];W) D C^([0, T*]; 7^*-^), we 
need the following Lemma. 

Lemma 2.11 (Continuity). AssumeU is a solution to (2.9) with T G W, and 
assume {U,dtU) satisfies {U,dtU) G C^„([0,T];H" x H'-^) n L~([0, T]; x 
H'-^), s > n/2 + 1. Then {U,dtU) G C([0,T];H" x H'-^). 

Proof. First recall that continuity on [0, T] is equivalent to right and left con- 
tinuity at each t G [0,T]. Changing the direction of time gives an equation of 
the same type, so it is sufficient to prove strong right continuity at t G [0,T]. 
By a reparametrization, there is no loss of generality in assuming t = 0. 

We know from Lemma 2.10 that t ^ {U,dtU) is weakly continuous. In 
order to prove strong continuity, we will use the fact that if Wm — ^ weakly, 
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then w„i w strongly if ^ II'^^IIj cf. [14, §12, exercise 3]. By fixing 

{g , N , X)p( (t)] in the definition of £, we may define a norm 

\\\U'\\\l., = £{t;{Dr-'U') (2.31) 

which for fixed t, is an equivalent norm to || • By (2.18), (2.15) and Lemma 
2.6, 

\\\m\\\s;t^\\\m\\\s;0 SSt\0. 
By compact imbedding, U G C{[0,T];H^ ) for s' < s, and hence we have g G 
C{[0,T];H^'). Prom this fact it follows easily that 

ll|W(i)||U^|||W(0)|||,;o ast\0. (2.32) 

and hence U is right continuous in Ti.^ at t = 0. The corresponding property 
for dtU follows from equation (2.9). □ 

By Lemma 2.11, we have U G C([0, T*]; H") n C\[0,n];n'~^). It fol- 
lows from (2.9) and the assumptions on the map U {g, N, X, J^) that U G 
CJp+\H'), where m is the order of regularity of (2.9). The contraction property 
used in the proof of Lemma 2.9 shows that U is the unique solution to (2.9). 

2.3. Cauchy stability. In this section we will give a proof of Cauchy stability 
following [2, 3]. We will give the proof only for the case s = k, k integer, 
k > n/2 + 1, and assuming that (2.9) is regular of order m = k — 1. The proof 
is easily adapted to noninteger s > n/2 + 1 and general m. 
Let II • 11^ denote || • ||-^« for £ integer. Introduce the norm 

\mm\k= E 

0<j<fe-l 

Define the spaces 

CtCH') = no<j<k-iC^i[0,T];n''-^), 
C^riv!') = r\o<j<k-iW^'\%T]-v!'-^), 

with norms 

|||Z^|||fc,T= E WWci {[0,T\-H''-i) 

0<j<k-l 

{{U))k,T = E ll^llw^JM([0,Tl;Wfe-J) 
Q<j<k-1 

Consider the linear problem 

L[h\U = T, U\^^^=U^ 

where we use h to denote the coefficients g,N,X of L. We will use the con- 
vention that |||/i|||^,T and {{h))^^ denotes the norm defined analogously to the 
above but using i?*^ instead of 7^^. From the energy estimate, Lemma 2.7 we 
get for 2 < ^ < A;, ^ integer, 

iiiwiii,,T < cm\\,,T)m% + iii.F(o)iii,_i + mu,T) (2.33) 

Using the density of C°° in Sobolev spaces, the proof of the following approxi- 
mation Lemma straightforward. 
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Lemma 2.12 (Approximation). Given any W e H^.T e C^iU^), i > 2, £ 
integer, there are for any e > 0, W° G C°°,J^e £ C°° such that 

||W°-W°||,<e 
111^(0) -^,(0)|||,_i< 6 
((^-^e)kT<e 

□ 

2.3.1. Perturbation estimate for i < k — 1. Consider the linear problems 

L[h]U = T, U\^^^=U^ 
L[h'\U' = F', U'\^^^=U'^ 

In the following write L = L[h], L' = L[h'] and assume h,h' G Ct{'H''). 
In applying the energy estimates in the following we will let C be a constant 
depending on T as well as on A[^], |||/i|||fc,Tj |||/i'|||fc,T- All of these quantities 
are under our control, using the apriori estimates. 

We calculate 

L'{U' - Ue) =J^'-Te + {L- L')Ue 

Let 

^^{e) = \\U% + mm\V-l + {{:Fe))^,T 
so that by (2.33), for l<k, 

\\\U,\\\t^T<C^t{e). 

Then, for £ < A: - 1, 

(((L - L')Ue))i,T < C\\\h - h'\\\k-i,T^e+i{e) 
Now we get from (2.33), for ^ < A; - 1, 

\\\u'-uMi,T<c{\\u'^-u% + \\\j''{o)-:FM\\\i-i 

+{{^' - .F,)),,T + \\\h- h'\\\k-i,T^i+i{e)] 

It is important to note that this works only for (. <k — 1. 

The energy estimate (2.33), gives when applied to W — We for £ < /c — 1, 

\\\U-U,\\\^^T<eC. 

Putting this together gives for £ < k — 1, 

\\\U-U'\\ V,T <C{e+\ \U'' -U% + \\ \m - J^'{Q) 1 1 

+{{J^-:F')),,T + \\\h-h'\\\k-i,T^i+i{e)] 
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2.3.2. Perturbation estimate for £ = k . Let 

B = dt-X 

so that 

LU = BU- JU 

where 

Let 5 = BU. Then 5 solves 
where 

= [S, J]U + Sj?^ 
We have the corresponding primed identities. 

Remark 2.3. For general s we have T under control in Tl^~^ , and this is also 
true forU^ . This observation allows one to generalize the proof to general s,m. 

From (2.34) apphed to 6 — S' , wc get 

\\\5 - 6'\\\k-i,T <C{e + \\U''-U%.i + (0) - .^'(0)|||fc_2 

+{{:F-^'))k-i,T + \\\h-h'\\\k-i,T^k{e)} 

The terms \\U'^ -U°\\k-i, |||J^(0)-/''(0)|||fc_2, {{T-P))k-i,T can be estimated 
in terms oi\\U'^ -U^\\k, |||.F(0) - .F'(0)|||fe_i, {{J^ - p))^,^. 
This gives 

\\\5 - 5'\\\k-i,T < C {e+ - U% + \mO) - J^'mWk-i 

+{{^ - :F'))k,T + \\\h- h'\\\k-i,T^k{e)} 

where now C also depends on ||Z^°||fc, 

We now use an elliptic estimate for W — W'. A computation shows 

J{U - U') = 5 - 5' + J'' - + {f - J)U' (2.35) 

Prom the definition of J and standard elliptic theory we get the estimate 

||w||fe<c(||Jw||fe_i + ||w|U-i) 

Hence (2.35) implies the estimate 

\\U-U'\\k<C{\\5-5'\\k-i + \\J'' -Hk-i + \\h-h'\\k-i\\U'\\k + \\U-U'\^^^^ 
Now calculate 

dt{U - W) =BU- B'W + XU- X'W 

= 5-5' + X{U - U') + {X- X')U' 
Iterating this estimate gives together with the above 

1 1 |W - W'll k,T < C (e + I |W° - I Ife + 1 1 \m - ^'(0) 1 1 U_i 
+{{^-J"))k,T + \\\h-h'\\\k-i,T^k{.e)] 
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2.3.3. Application to the nonlinear system. Consider the Cauchy problems 

L[U]U = J^[U], Ul^^=U°, 

We wish to estimate \\\U — Z^'|||fe,T in terms of W° — By the assumptions, 
if U^U' are close, there is a constant CLip (not to be confused with C^), such 
that 

\mum - nu'm\\\k-i < CLipiiw°-w'°iife-i, 

{{T[U] - J'p(']))k,T < CupiiU - U'))k,T, 

|||/i-/i'||U-i,r<CLip|||Z^-Z^'|||fe-i,T. 
By (2.36) we now get an estimate of the form 
\\\U-U'\\\k,T<C{e+\\U''-U'% 

+ {{U-U'))k,T +\\\U- U'\\\k-l,T^k{^)} . 

The term {(JA — U'))k,T can be eliminated from the right hand side by an ap- 
plication of the Gronwall inequality. This gives 

\\\U- U'\\\k,T < C (e + I |W° - U + 1 1 |W - W'l I \k-i,T^k{.e) } . (2.37) 

Now consider a sequence {U^} C H^, such that W° — W° in as a — oo. Let 
Ua solve 

By compact imbedding we can choose a subsequence so that 1 1 — Z^|||jfc-i,r ~^ 
0. It follows from (2.37), using the fact that e was arbitrary, that 

— ^ 0, as a — > OO. 

This completes the proof of Cauchy stability and point 1 of Theorem 2.3 is 
proved. 

2.4. Continuation. It remains to prove the continuation principle, point 2 
of Theorem 2.3. Suppose for a contradiction, (2.11) does not hold and that 
r+ < oo is the maximal time of existence. Then by the proof of local existence, 
there is a uniform lower bound for the time of existence with initial data U{t), 
t < T+. This contradicts T+ < oo. The fact that T+ is continuous follows 

from Cauchy stability. This completes the proof of Theorem 2.3. 

3. The Cauchy problem for the modified Einstein evolution 

equations 

Let n>2 and fixs>n/2 + l. The vacuum Einstein equations have special 
structure in dimension 2 + 1 which we do not make use of, but we allow n = 2 
here and in the following sections for completeness. 

Let P be the operator defined by 

PY' = AX' + RKxf - jCxV - 2V"*X"eXV„e„ - V^e„), 
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where V is given by (1.5). Let the operators B,E be defined by 
Bf = -Af+\kff, 

Ef = -2v^fki + vv^r + vk^^'iyLn - rD- 

With A given by 

(3.1) 

the defining equation (1.11) can be written in the form 

4^)-ri). (3.2) 



X ) 

The operator A is second order elhptic. 

Define V to be the set of symmetric covariant tensors (5, k) e x H^~^ 
such that 

5 is a Riemann metric (3.3a) 

The operators B, P : , are isomorphisms at {g, k). (3.3b) 

Suppose V is nonempty. For {g, k) € V, let Cnx = Cnx{9, k) be a constant so 
that for functions u and vectorfields Y, 

\\u\\l2 < Cnx\\Bu\\l2, (3.4a) 

\\Y\\l2 <Cnx\\PY\\l2. (3.4b) 

Let {g^,k^) G V be given. In this section we prove that the Cauchy problem 
for the system (1.10,1.4,1.6), with initial data 

(5,fc)[^^ = (<7°,fc°), (3.5) 

is strongly locally well-posed. We will refer to this problem as the CMCSH 
Cauchy problem with data (5°,^;°). 

Theorem 3.1. The CMCSH Cauchy problem with initial data {g^,k^) e V is 
strongly locally well-posed in C'^{W), k = [s\. In particular, there is a time of 
existence > so that the solution map {g^,k^) {g, k, N, X) is continuous 

Here can he chosen to depend only on C]\fx{g^,k^), ^[9^], Wo^Wh'' cLn-d 
Wk^W^s-i- In particular, T* can be chosen so that it depends continuously on 
(5°,A;°) eH^ X H^-^. 

LetTj^ he the maximal time of existence of the solution to the CMCSH Cauchy 
problem with data {g^,k^). Let Cnx = CNxUdit), k{t)) be defined by (3.4)- 
Then either r+ = 00 or 

limsupmax(A[3], ||L>3||z,oo, ||A:||z,oo,CArx) = 00 

To prove Theorem 3.1 we will show that the CMCSH Cauchy problem is 
quasi-linear hyperbolic in V, regular of order [s\ — 1. The result then follows 
from Theorem 2.3. The following Lemma gives the basic estimates used in the 
proof of Theorem 3.1. 
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Lemma 3.2. The set V defined by (3.3) is an open subset of x Let 
N,X be the solution to (1-11) w.r.t. (g, k) and let g be defined by (1.1) in terms 
of {g,N,X). Let the map a : {g,k) i-^ {N,X) be defined by solving (l-ll)- 

There is a constant Cl > depending only on Cjvx, ||D^||ioo such that 
the following holds. 

1. S^/^^(5°,/c°) C V, where B'y^^{g^ ,k'^) is the ball in W of radius 1/Cl 

centered at {g^,k^). 

2. The Frechet derivative Da satisfies 

\\Da{g,k){g',k')\\H^<CL\\{g',k')\\n^-i 

for {g,k)eBl^^^{g\k^). 

3. a has the Lipschitz property 

Mg,k)-a{g',k')\\Hr+i < CLUg,k) - ig',k')\\n^. 

for all (g, k), {g\ k') G k^), l<r<s. 

4. For integers m, j, ii, such that m = [sj — 1, 1 < j < m, 1 < £i, J2i — 
the Frechet derivative D^a is a Lipschitz map from B{j^^{g^,k^) to the 
space of multilinear maps 

3 
1=1 

Remark 3.1. Note that the Lemma gives an estimate for N,X e H^'^^. This 
extra regularity for N, X will be important later in the proof of Theorem 3.1. 

Proof. Consider the second order elliptic system for {N,X) given by (1.11). 
In this system the zeroth order coefficients depend on the Ricci tensor R--' and 
derivatives of the vector field given by (1.5), and thus, at first glance, appear 
to contain second derivatives of gij, i.e. terms in H^~^. In order to prove that 
{N, X) G i?*"^^ X H^'^'^ for {g, k) G V, we must prove that these second derivative 
terms cancel. 

It is clear from (1.11) that it is sufficient to restrict our attention to the 
system 

^ AY' + R'fYf - CyV' 
where is given by (1.5). We have 

AFe + RjYf = g'^VaiVbYc + VeYb - VmY'^gbc) 

We will compute in local coordinates. Writing VfoFc + VftFc = ^Y9bc in terms 
of di gives 

^ygbc = Y^d^gbc + QbfdcY^ + gcfdbYf. 
Using this together with VmY"^ = \g™'^CYgmn shows after some manipulations 

that the terms in AYc+RjYf containing second derivatives of gab can be written 
in the form 

^g'^^'Y^dm {da9bc + dbgac - dcgab) . (3.6) 
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On the other hand, using the exphcit form of the Christoffel symbol, 

= \9'\dj9ki + dkQji - diQkj) 

and the form of in local coordinates (1.5b) one sees that (3.6) exactly cancels 
the terms in gcdJO-yV^ containing second order derivatives in gij. 

Let the operators B,P,Aheas in (3.1). It follows from the fact that A is lower 
triangular that if B,P : — > are isomorphisms, then also A : — > is 
an isomorphism. 

With u = {N,X), (3.2) and hence (1.11) is of the form 

A{g,Dg,k)u' =F\g,Dg,k) (3.7) 

where F is a smooth function of its arguments and ^ is a second order elliptic 
system of the form 

A{g, Dg, k)u' = + b'^g, Dg, k)diu"' + c^g, Dg, k)u"', (3.8) 

where b, c arc smooth functions of their arguments. It now follows from standard 
elliptic theory, cf. [17] that an inequality of the form 

\\u\\Hr+2 <C{\\Au\\Hr + \\u\\Hr), 0<r<s-l, (3.9) 

holds. Uniqueness together with compact imbedding implies that the lower 
order term can be eliminated from this inequality. We sketch the standard 
argument for this. 

Suppose that an inequality of the form ||n||j^r+2 < C||An||j/r does not hold. 
Then there is a sequence {ui}^^ with ||tii||^r+2 = 1, and — > as 
i ^ oo. Since M is compact, is compactly imbedded in H^, and therefore 

there is a subsequence {uj} which converges to some which by construction 
satisfies Au^ = 0. By (3.9), in fact G H''^'^. The existence of a solution 
to Au = contradicts the assumption that equation (3.7) has unique solutions. 
Thus, an inequality of the form 

lh||ijr-+2 < Cp-ullifr-, 0<r<s-l, (3.10) 

holds. It follows from the commutator estimates as in Lemma 2.7, that the 
constant C in (3.10) depends only on C^. 

It is convenient to use the notation h = (g,k), = x H^~^. We 
write A[h] = A[g,Dg,k], F[h] = F{g,Dg,k) etc. To prove that C depends 
continuously on h eW, let h' be close to h in W and calculate 

ll^I/i'lnlli/-- > PWulli/r - \\{A[h']-A[h])u\\Hr 

> C~''"||u||^r+2 — Cl||/l' — ^||?is||'u||^r+2 

which by choosing h' sufficiently close to h shows that the estimate |]';i||/^r+2 < 
(7'||Au||/fr holds in a neighborhood of g with a uniform constant C arbitrarily 
close to C. In particular, for h' close to h, A[h'] is an isomorphism. From this 
follows easily the first part of the Lemma and point 1. 

Point 2 follows from elliptic estimates for (3.7). The Lipschitz property, 
point 3 for the map {g, k) {N, X) is a consequence of 2 and the mean value 
inequality. The proof of point 4 is straightforward, the basic tool being the 
following product estimate which is a consequence of the product estimate II. 
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Let r > n/2, ii > 1, i = 1, . . . ,j, = m, m < r. Then multiplication is a 

bounded multilinear map 



H''-^^ 



□ 



Proof of Theorem 3.1. We will show that Theorem 2.3 applies to the modified 
Einstein evolution equations (1.10,1.11). In order to do this we must first write 
the system (1.10,1.11) with initial data {g,k)\^__^^^ = {g^,k^), as a quasi-linear 
hyperbolic system (in the sense of Definition 2.2) of the form 

Let 

= 9i3^ (3.11a) 
Vij = -2kij, (3.11b) 

Then Uij,Vij are symmetric 2-tensors. We expand the Lie derivatives CxUij, 
^xVij in terms of V, 

and similarly for CxVij. Similarly ViViN can be written in terms of V as 

ViVjN = ViVjN - ^g^'^iVigjm + Vjgim - ^mgji)ViN. 

The Ricci tensor Rij of gji is quasilinear elliptic, up to a gauge term. This 
is seen from the identity 

Rij = ——Aggij + Sij + Sij 
where Sij[g, dg] is given by 



Sij = 7ii9li9"'"'R^mjn + 9lj9"^'^R^min) 



-j{9li9""' mjr. 
+ ^g^'g^' S ^jgns^Wim + ^igml^sgjn - ^^jgns^igml 



+^m9il^ngjs — Vm5i«Vsffjn| 



l^Vigij, 



where 



* y \^ rnn mnJ 

In the computation giving Sij, is it convenient to make use of the fact that Sij 
is a tensor, and work in a local coordinate system with V'^-j^ = at the center 
of coordinates. 
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Let Sij[u,du] be the expression corresponding to Sij[g,dg]. Now we may 
write the system that shows up in the proof of local existence in the form 



with 



N -25,;. 



dtUij = Nvij + CxUij 

dtVij = NAgUij + CxVij + Fij 



' ran , mn 



(3.12a) 
(3.12b) 



+ 2 



VmUji)VlN 



(3.13) 



where ■u"*" = {u ^)mn- 

In order to apply the energy estimates as presented in section 2, we further 
expand the Lie derivative in terms of V. This gives the system in the form 

dtUij = Nvij + VxUij + J='iij (3.14a) 
dtVij = NAgUij + Vxvji + J^2ij (3.14b) 

^2ij = Fij + vijViX^ + ViiVjX^ 

with Fij given by (3.13). Finally, assuming that the defining equations for N, X 
have unique solutions, let N = N[u,v\, X = X[u,v\ be given by (1.11). Then 
(1.10) takes the form 



where 



dfU = Nv + Vxu + ^1 [u, v] 
dfV = NAgU + Vxv + J^2 [u, v] 



(3.15) 



Introduce the notation U = {u,v), then g = g[U],k = k[U] and 
{J^i[U],J^2[U]). Further, via (1.11) we have N = N[U],X = X[U]. 
Then we can write equation (3.15) in the form 



L[U]U = J^[U], U\ 



t=to 



(3.16) 



with L given by (2.8), andU^ given in terms of {g^, k^) by (3.11). Clearly (3.16) 
is a system of the form considered in Definition 2.2. From the definition of 
and Lemma 3.2, it follows that the CMCSH Cauchy problem with data in V is 
quasilinear hyperbolic, regular of order m = \_s\ — 1, and hence Theorem 3.1 
follows from Theorem 2.3. □ 



Remark 3.2. The proof that T : — > H.^ is the only place in the proof of 
Theorem 3.1 where the "additional regularity" N,X e x is used, of 
Remark 3.1. 
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4. Evolution of gauges and constraints 

Let n > 2 and fixs>n/2 + l. Let ^ be a fixed C°° metric on M with Levi- 
Civita covariant derivative V. Introduce the constraint and gauge quantities 

{A,F,V\D') by 

A = tYk- t, (4.1a) 
y*^ = 5*Je'=(V,e, -V.e,) (4.1b) 

F = R+{trkf -m"^ -ViV\ (4.1c) 
A = VitrA: - 2V'^kmi- (4. Id) 

Let also 

<5«.- = ^(V^F,+V,y,), (4.2) 

A calculation shows the constraint and gauge quantities (^, -F, F*, -D*) satisfy 
a hyperbolic system when (gi, k, N, X) solve the modified evolution equation 
(1.10-1.11). Using an energy estimate for this hyperbolic system, we will prove 
that if {A,F,V\D'^) = for the initial data ig^,k^), then (A,F,V\D'') = 
along the solution curve {g, k, N, X) with 

(g, k) e no<i<rn < C^((ro, Ti); H'-^ x H'-^'^), (4.3a) 
(TV, X) G no<^<™C^((ro, Ti); H'+'-' x H'+'-'), (4.3b) 

where m = [s\ — 1, constructed in Theorem 3.1. In computing the time deriva- 
tives of the constraint quantities A, F, V, D, we note the fact that the effect of 
Lie dragging gij, hij by X, is that the quantities A, F, V, D are also Lie dragged. 
We have using (1.10-1.11) and (4.1), 

dtV' = CxV' + ND^ (4.4a) 

dtA = CxA + NF (4.4b) 

dtF = CxF + ViND' + 2NFtvk + 2Nk'^dij + V^VfiNtrk) 

+ NAA + 2ViNV'A (4.4c) 

dtDi = CxDi + ViNF + NtikDi + N{AVi + RifV^) 

+ NtikViA + 2V'^N6mi (4.4d) 

In doing these computations we have used the expressions for the Prechet deriva- 
tives of of r and R, 

DT%.h = ^g'^'iVjhkm + ^khjm - ^mhjk) (4.5) 

DR.h = -A/i> + V'V^hij - Rijh'^ (4.6) 

and the identities 

-2AVjX^ + V V^' {ViXj + VjXi) = VmRX"" + 2R'"'ViXm (4.7) 



(4.8) 
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The leading order terms in (4.4) are 

dtA ^ NF 
dtF ^ NAA 
dtV' ^ ND' 

dtDi^N{AVi + RifVf) 

Using the product and composition estimates stated in section 2 and the defi- 
nition of Rij in terms of gij, one finds that at a Riemann metric gij, the map 
gij I— Rj^j is smooth and satisfies 

\\Rij\\Hs-2<C{A[g])\\g\\Hsil + \\g\\H^), 

for s > n/2 + 1. 

Define the energy 6 = 6i + 62 hy 

£1 = 11 {\A\' + \VAf + \Ff)^,, 

and 

£2 = ^ l{\V\^ + \VVf + \D\^)^ig 

The foUowing Lemma gives the energy estimate required for proving that the 
Einstein vacuum constraints and the CMCSH gauge is conserved by the modi- 
fied Einstein evolution equations. 

Lemma 4.1. Let {g,k,N,X) be a solution to (1.10,1.11), satisfying the regu- 
larity condition (4-3). There is a constant C = C{g,k,N) so that 

\dt£\ < CS 

Proof. We compute dt£ using (4.4). Due to covariance, the Cx terms in (4.4) 
can be dropped. From the assumptions, it follows that dtgij G H^~^ C L°°. 
Therefore, we only need to consider the terms in dfS involving dfA, dtF, dtV^, 
and dtDi. Using (4.4) and performing a partial integration, it is easy to check 
that 

\dt£i\ < CS. 

It remains to consider dtS2. It is straightforward to show dt J \V\'^iJ.g < C£. It 
remains to consider 

dtl J iViV^VkV'g"'gji + DiDjg'^)fig 

The only terms which need to be considered in detail are those where dt hits 
r^^, or Di. The term involving dfT-'- is of the form 

tilt i/fit 

j {^tTlJV^V,V^g^>^g,l^ig (4.9) 
The terms involving dtV^ , dtDi yield, after a partial integration, the expression 

j NRifV^D^Hg (4.10) 

Prom the assumptions, we have dtTj^ and Rif bounded in H^~'^. Further, V 
and D are bounded in and L^, respectively, by £. The expressions (4.9) 
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and (4.10) are of the form J uvw with u € H^^'^, v G H^, w E L^. By the 
product rule (2.6), with ti = s — 2, t2 = 1, p = 2, we have since s > n/2 + 1, 
||it^^||L2 < C'||ti||j:^s-2||w||j:^i. Au apphcatiou of the Cauchy inequahty gives the 
estimate 

I J UVWflgl < C\\u\\jjs-2\\v\\hi\\w\\i2. 

Together with the above this gives \dt£2\ < C£. This completes the proof of 
Lemma 4.1. □ 

Lemma 4.1 and the Gronwall inequality now shows that if £^ = initially, then 
6 = along the solution curve. We now have 

Theorem 4.2. Let {g^,k^) € V, and assume {M,g^,k^) satisfy the Einstein 
vacuum constraint equations (1-4) cind the gauge conditions (1-6), i.e. 

{A,F,V,D)=0. 

Then the space-time metric g defined in terms of the solution {g, k, N, X) of the 
CMCSH Cauchy problem, is a solution of the Einstein vacuum equations. □ 

5. Isomorphism property 

Let n > 2 and fix s > n/2+ 1. In this section, we use local coordinates unless 
otherwise stated. Define Q to be the set of {g, k) G x H^~^ such that 

g \s a Ricmann metric (5.1a) 
(A, F, L>) = at (5, A;) (5.1b) 

Thus {g,k) G Q precisely when {g,k) satisfies the constraint equations (1.4) 
and the gauge conditions (1.6). In this section we will work in the CMC time 
t = tigk. For initial data {g^, k^) with tg = tigok^, let (T_,T+), < tg < T+ 
be a maximal existence interval, defined by analogy with the notion of maximal 
existence time. 

Theorem 5.1. Assume that {g^,k^) G Q and that g has negative sectional 
curvature. Then the CMCSH Cauchy problem is strongly locally well-posed 

in C^{J-i^), k = [sj and the Lorentz metric g constructed from the solution 
{g, N, X) is a vacuum solution of the Einstein equations. Further, the following 
continuation principle holds. Let to = tigok^ < and let {T-,T^), T_ < io < 
be a maximal existence interval for the CMCSH Cauchy problem in CMC 
time t = tigk. Then either (T_,r+) = (— oo,0) or 

limsup(A[5] + ||i:>5||L- + H^Hl-) = 00 
ast /Tj^ or ast\ T_. 

We begin by considering the operator P defined by 

PY' = AY' + R'fY^ - CyV' - 2V"'F"(r^„ - f j„„) 

We will need some material concerning vector and tensor fields along maps. 
Consider a map (f) : {M,gij) —>■ {N,hap). We use latin indices for coordinates 
on M and greek indices for coordinates on N. The bundles of tensor fields 
along (f), ®^T*M (g) (f)~^TN have natural connections which we denote by D. 
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Here (p^'^TN is the pullback of TN to M along (j). We work this out in local 
coordinates for = 0, 1, 2. 

Div^ = div^ + "'rl^v'^dicj)^ (5.2) 

ACJ = - ""^ijCk + ^r^l^ej'S,/ (5.3) 
AC], = d,q, - ^vl,Q, - ^Ti,Q + ^rl^^^d,^ (5.4) 

where if the left hand side is evaluated at a; G M, all objects on N are evaluated 

at (p{x). 

On the bundles ®^T* M®(t)-'^TN , there is a natural inner product (•, •) which 
is local coordinates is given by 

{v,w) = v°'w'^hag 

and similarly for higher order tensors. It is straightforward to check that D is 
metric w.r.t. (•, •). Further, letting Ajjv"" = g^^DiDjv"' the identity 

/ {ADV,w)iJ,g = - {Dv,Dw)ng 
Jm Jm 

holds, and thus the operator Ad is self-adjoint with respect to the pairing 

{V,w)flg 



/( 

JM 



Lemma 5.2. Assume g has negative sectional curvature and V = 0. Then 

P : H^~^^ H^~^ is an isomorphism. 

Proof. Let : M — > M be a diffeomorphism. Then 

This follows from that fact that the difference of Christoffel symbols transforms 
as a tensor, or by a direct computation using the identities 

Let y be a vector field on M and let (l)s defined hy ds4>s = Y ° 4>s, '/'o = Id, be 
the flow of y. A computation shows 

and similarly for ^'^*"^^T\-. Recall that 



Now consider the identity (5.5) with (p replaced by ^g- Differentiating with 
respect to s and evaluating at s = 0, gives the identity 

pyi ^ ^^n^^Jj^yi ^ Ri^.^Y^) (5.6) 
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where the indices on the Ricmann tensor R oi g are raised and lowered with g. 
Using the definition of V in terms of T we have 

5™"V^Vnr^ = g"'''{dm^nY' - t^rnn^rY' + t^J/nY^) 

= g'^'^idJlnY' - Tl^JlrY' + tiJ/r^Y^) + V'V rY' 

where we used the definition of V to replace F by F. 

Now we will think of y as a vector field along the map Id : (M, g) (M, g). 
For clarity we will use lower case greek indices for objects associated with g. 
Recalling the definition of the connection D and the operator Ad, we see that 

5"'"V„,V„y° = Ai^y^^ + VDrY'^ 

Hence equation (5.6) takes the form 

Py° = Aoy" + g'^^'R^p^Y^ + V^DrY"" 

Now assume as in the statement of the Lemma that V = 0. Then 

[{PY,Y)fig = - [ {DY,DY)iig+ f g'^''Ram/3nY''Y(^i^g (5.7) 

J JM JM 

The quadratic form (Y, W) g"^"" Rami3nY°'W'^ is symmetric and when g has 
negative sectional curvatures then there is a A > so that 

g'^'^'Ran^pnY'^YP < -X^\Yf 

It follows that an estimate of the form 

\\Y\\L2<C{A[g])\\PY\\L. 

holds, and the Lemma now follows from the fact that P is second order elliptic 
of the form (3.8), cf. the proof of Lemma 3.2. □ 

Proof of Theorem 5.1. Let P be as above, let the operators B,E,A be defined 
as in section 3. By the proof of Lemma 3.2, A is an elliptic second order operator 
of the form (3.8). Therefore, in order to prove that A is an isomorphism, it 
is enough to prove that ker A = ker^* = {0}. It follows from the maximum 
principle that B is an isomorphism as long as tvgk ^ 0. By Lemma 5.2, P is an 
isomorphism. In view of the fact that A is lower triangular, the isomorphism 
property now follows from the isomorphism property for B and P. 

We have now proved that Q G V and hence the conclusion of Theorems 
3.1 and 4.2 apply. It remains to prove the continuation principle. In view of 
Theorem 3.1, it is enough to estimate Cnx in terms of A[^] and (tr^A;)""^. By 
(5.7) it follows that ||y||L2 < C{A{g))\\PY\\L2. Similarly, we have 

{BfJh^> f \Vff + \k\'f>C{A{g))^^^\\f\\l, 

JM <J 

which gives the corresponding estimate for B. Now we have bounded Cjvx 
in terms of A[^] and (tr^A;)"^. Finally, to see that T+ < 0, note that as M 
admits a metric g with negative sectional curvature, it admits no metric with 
nonnegative scalar curvature, cf. [7, Corollary A, p. 94]. In case n = 2, this 
is a consequence of the uniformization theorem. Now suppose that T+ > 0. 
Then tr^A; = must occur during the evolution, at which instant the scalar 
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curvature of g is nonnegative by the Hamiltonian constraint (1.4a), which gives 
a contradiction. Therefore we have < 0. This completes the proof of 
Theorem 5.1. □ 
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